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Abstract. We review the present status of the nucleon effective mass splitting puzzle in asymmetric matter, with contro- 
versial predictions within both non-relativistic and relativistic approaches to the effective in medium interactions. Based on 
microscopic transport simulations we suggest some rather sensitive observables in collisions of asymmetric (unstable) ions 
at intermediate (RIA) energies: i) Energy systematics of Lane Potentials; ii) Isospin content of fast emitted nucleons; iii) 
Differential Collective Flows. Similar measurements for light isobars (like 3 // He) could be also important. 



ISOSPIN MOMENTUM DEPENDENCE IN SKYRME FORCES 



We start from a discussion on the isospin dependence of these widely used non-relativistic effective interactions 1 1]. 
In a Skyrme-like parametrization the symmetry term has the form: 
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with the function C(p), in the potential part, given by: 
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with a > and the usual Skyrme parameters. We remark that the second term is related to isospin effects on the 
momentum dependence 01 • 

In the Fig.^left) we show the density dependence of the potential symmetry term of various Skyrme interactions, 
SIIL SGII, SKM* (see 0| and refs. therein) and the more recent Skyrme-Lyon forms, SLya and Slyb (or SlyA), see 
We also separately present the local and non-local contributions, first and second term of the Eq.l|2j. We clearly 
see a sharp change from the earlier Skyrme forces to the Lyon parametrizations, with almost an inversion of the signs 
of the two contributions, |5]. The important repulsive non-local part of the Lyon forces leads to a completely different 
behavior of the neutron matter EOS, of great relevance for the neutron star properties. Actually this substantially 
modified parametrization was mainly motivated by a very unpleasant feature in the spin channel of the earlier Skyrme 
forces, the collapse of polarized neutron matter, see discussion in |8|, |2j, |4j, |gfl . In correspondence the predictions on 
the isospin effects on the momentum dependence of the symmetry term are quite different, see Fig^(left). A very 
important consequence for the reaction dynamics is the expected inversion of the sign of the n/p effective mass 
splitting. 



Effective masses in neutron-rich matter 



In asymmetric matter we consistently have a splitting of the neutron/proton effective masses given by: 
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FIGURE 1. Left: Density dependence of the potential symmetry term for various Skyrme effective forces, see text. The bottom 
right panel shows the kinetic contribution. Right: Density dependence of the neutron/proton effective mass splitting for various 
Skyrme effective forces, see text. The asymmetry is fixed at I = 0.2, not very exotic. 
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The g\ , g2 coefficients are simply related to the momentum dependent part of the Skyrme forces: 

Si = -^2[fi(2+xi)+f 2 (2+x 2 )] 
An 

g2 = -U[f2(l+2x 2 )-?i(l+2*i)] (4) 
An 

This result derives from a general q — structure of the momentum dependent part of the Skyrme mean field 

U qM D = m(gip+g 2 p q )E (5) 
where E is the nucleon kinetic energy, while the total field seen by the g-nucleon has the form |7] 

TTf nn tA tj gift 2 T g 2 h 2 x q de Loc (Pot) 

U q {p,p q ,k)=U qMD + —^—+ 2 H + — (6) 

where h 2 T is the kinetic energy density and eu )c (Pot) the local part of the potential energy density. 

In the Fig. fright) we show the density behavior of m* p in neutron rich matter I = 0.2 for the same effective 
interactions. From the Eqs.(|3]0} we see that the sign of the g 2 univocally assigns the sign of the splitting, i.e. g 2 < 
gives larger neutron masses m* > m* p while we have the opposite for g 2 > 0. 

In the Table [l\ we report some results obtained with various Skyrme forces for quantities of interest, around 
saturation, for the present discussion. We show also the E — slope of the corresponding Lane Potential, see later, 
simply related to the isospin dependent part of Eq.(0. For the effective mass parameters of Eq.(|3} we observe that 
while the g\ coefficients are always positive, corresponding to a decrease of the nucleon mass in the medium, the 
isospin dependent part shows different signs. In particular we see that in the Lyon forces the g 2 values are positive, 
with neutron effective masses below the proton ones for n-rich matter as shown in Fig.^(right). 



TABLE 1. Properties at saturation 



Force 


Sill 


SGII 


SkM* 


SLya 


SLyA 


SLyl 


g x (10- 3 )(MeV-'/m 3 ) 


+3.85 


+3.31 


+3.53 


+ 10~ 5 


+ 1.67 


+ 1.70 


g 2 (\0- 3 )(MeV- l fm 3 ) 


-3.14 


-2.96 


-3.50 


+5.78 


+2.53 


+2.76 


PoO* -3 ) 


0.150 


0.1595 


0.1603 


0.160 


0.1595 


0.1581 


a 4 (MeV) 


28.16 


26.83 


30.03 


31.97 


32.01 


32.01 


C(p )(MeV) 


31.72 


29.06 


35.46 


39.40 


39.42 


39.42 


E — slope(LanePot .) 


-0.22 


-0.21 


-0.26 


+0.43 


+0.19 


+0.20 



We note that the same is predicted from microscopic relativistic Dirac-Brueckner calculations and in 

general from the introduction of scalar isovector virtual mesons in RMF approaches fl^. ll3ll. At variance, non- 
relativistic Brueckner-Hartree-Fock calculations are leading to opposite conclusions 1141 llfll . We remind that a 
comparison between relativistic effective (Dime) masses and non-relativistic effective masses requires some attention. 
This point will be carefully discussed later. 



Energy dependence of the Lane Potential 



The sign of the splitting will directly affect the energy dependence of the Lane Potential, i.e. the difference between 
(n,p) optical potentials on charge asymmetric targets, normalized by the target asymmetry 0. From the Eqs. (I5I6> 
we obtain the explicit Skyrme form of the Lane Potential: 
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where C(po) gives the potential part of the 04 parameter in the mass formula. We see that the E — slope has just the 
sign of the g2 parameter, and so we have opposite predictions from the various Skyrme forces analysed here. The 
change in the energy slope is reported in the last row of the Tabled The difference in the energy dependence of the 
Lane Potential is quite dramatic. 

The second term of the Eq.Q is also interesting. In the case g2 > (positive slope) decreases the starting zero 
energy point, while in the case g2 < (negative slope) it represents a small correction to the symmetry energy C(po). 
We then expect to see a crossing of the two prescriptions at very low energies, i.e. low momentum nucleons will see 
exactly the same Lane potentials. 

An important physical consequence of the negative slopes is that the isospin effects on the optical potentials tend 
to disappear at energies just above 100 MeV (or even change the sign for "old" Skyrme-like forces). Unfortunately 
results derived from neutron/proton optical potentials at low energies are not conclusive, since the effects 

appear of the same order of the uncertainty on the determination of the local contribution. Moreover at low energies 
we expect the crossing discussed before therefore for a wide energy range we cannot see differences. More neutron 
data are needed at higher energies, in particular a systematics of the energy dependence. 

We can expect important effects on transport properties ( fast particle emission, collective flows) of the dense and 
asymmetric NM that will be reached in Radioactive Beam collisions at intermediate energies, i.e. in the RIA energy 
range. 



SYMMETRY ENERGY IN QUANTUM-HADRO-DYNAMICS 

The QHD effective field model represents a very successful attempt to describe, in a fully relativistic picture, 
equilibrium and dynamical properties of nuclear systems at the hadronic level Here we focus on the 

dynamical response and static (equilibrium) properties of Asymmetric Nuclear Matter (ANM) in a "minimal" effective 
meson field model, with two meson (scalar and vector) contributions in each isospin channel, (a, Co) for the isoscalar 
and (5,p) for the isovector part, see details in ||5|] and refs. therein. In particular we will discuss the effects of the 
isovector/scalar 8 meson, since this is directly related to the nucleon effective mass splitting of interest here. 



The symmetry energy in ANM is defined from the expansion of the energy per nucleon E(pb,I) in terms of the 
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In the Hartree case an explicit expression for the symmetry energy is easily derived that can be reduced to a very 
simple form, with transparent 5-meson effects I12lll3ll : 
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where k.F is the nucleon Fermi momentum corresponding to p B , Ep = J (hp + m* 2 ) and m* is the effective nucleon 
mass in symmetric NM, m* = m + 'L s (a), with the isoscalar, scalar self-energy E s (a) = —fpPs (p.v scalar density). 
fi = (^-) 2 , i = <?, (0, 8,p), are the effective meson coupling constants |24j]. 

We see that, when the 8 is included, the empirical 04 value actually corresponds to the combination \f p — fsijrj;) 2 ] 
of the (p , 8) coupling constants. Therefore if fg ^ we have to increase correspondingly the p -coupling. 

Now the symmetry energy at saturation density is actually built from the balance of scalar (attractive) and vector 
(repulsive) contributions, with the scalar channel becoming weaker with increasing baryon density. This is clearly 
shown in Fig|2jleft). This is indeed the isovector counterpart of the saturation mechanism occurring in the isoscalar 
channel for symmetric nuclear matter. 
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FIGURE 2. Left: p- (open circles) and 8- (crosses) contributions to the potential symmetry energy. The solid line is a linear 
extrapolation of the low density 8 contribution. Right: Total (kinetic + potential) symmetry energy as a function of the baryon 
density. Dashed line (NLp). Dotted line (NLpS). Solid line NLHF, only Fock correlations MSI . 

In FigOright) we show the total symmetry energy for the different models. At subnuclear densities, p B < po, in 
both cases, NLp and NLp 8, from Ea.l llOt we have an almost linear dependence of E sym on the baryon density, since 
m* ~ Ef as a good approximation . Around and above po we see a steeper increase in the (p + 8) case since m* /Ef is 
decreasing. 

In conclusion when the 5-channel is included the behaviour of the symmetry energy is stiffer at high baryon density 
from the relativistic mechanism discussed before. This is in fact due to a larger contribution from the p relative to the 
8 meson. We expect to see these effects more clearly in the relativistic reaction dynamics at intermediate energies, 
where higher densities are reached. 



Nucleon Effective Mass Splitting 



An important qualitatively new result of the 5-meson coupling is the n/p-effective mass splitting in asymmetric 
matter, |5] and refs. therein, Dime mass, see later: 

m* D {q)=m+L s (a)±fsPs3 , (11) 

where + is for neutrons. Since P53 = p$ p — Psn, in n-rich systems we have a neutron Dime effective mass always 
smaller than the proton one. We note again that a decreasing neutron effective mass in «-rich matter is a direct 
consequence of the relativistic mechanism for the symmetry energy, i.e. the balance of scalar (attractive) and vector 
(repulsive) contributions in the isovector channel. 



Dirac and Schrodinger Nucleon Effective Masses in Asymmetric Matter 

The prediction of a definite m* D (n) < m* D (p) effective mass splitting in RMF approaches, when a scalar 8 — like 
meson is included, is an important result that requires some further analysis, in particular relative to the controversial 
predictions of non-relativistic approaches. 

Here we are actually discussing the Dirac effective masses m* D (n,p), i.e. the effective mass of a nucleon in the 
in-medium Dirac equation with all the meson couplings. The relation to the Schrodinger effective masses m* s (n,p), 
i.e. the "k-mass" due to the momentum dependence of the mean field in the non-relativistic in-medium Schrodinger 
equation is not trivial, see i25ll2lll26H27il . We will extend the argument of the refs. I26ll27ll to the case of asymmetric 
matter. 

We start from the simpler symmetric case without self-interacting terms. The nucleon Dirac equation in the medium 
contains the scalar self-energy T. s = —foPs an d the vector self-energy (fourth component) = f a pB and thus the 
corresponding energy-momentum relation reads: 

(e + m - E ) 2 = p 2 + (m + E s ) 2 = p 2 + m* D 2 (12) 

i.e. a dispersion relation 

e = -m + Lo+Jp 2 +m* D 2 (13) 



From the total single particle energy E = E + m expressed in the form E — \A 2 + m 2 , where is the relativistic 
asymptotic momentum, using Ea. (ll2> we can get the relation 
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^ + L s + L Q + ^(T. 2 -I 2 ) + ^e = ^ + U eff (p B ,p s ,e) (14) 
Zm 2m m 2m 

i.e. a Schrtidinger-type equation with a momentum dependent mean field that with the dispersion relation Eq . (1 1 3i is 
written as 

Ueff = S, + J— + Eq) + — V(P 2 + m* 2 ) 
2m m 
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m 2m 2m m n 



The relation between Schrodinger and Dirac nucleon effective masses is then 
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Since at saturation the two self-energies are roughly compensating each other, E s + £0 — — 5QMeV the two effective 
masses are not much different, with the S — mass slightly larger than the D — mass. 



In the case of asymmetric matter, neutron-rich as always considered here, we can have two cases: 
• Only p meson coupling 

Now the scalar part is not modified, we have the same scalar self energies T, s for neutrons and protons and so the 
same Dirac masses. The vector self energies will show an isospin dependence with a new term T,fpPB3, repulsive 
for neutrons (— sign, since we use the definition Pjb,s)3 = P(B.s)p — P{B,s)n)- As a consequence we see a splitting 
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at the level of the Schrodinger masses since Ea. (ll6> becomes 

in 

mi(n,p)=m* n — — — (17) 
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in the direction of ms(n)* < ms(p)* (here and in the following upper signs are for neutrons). 

• p + 8 coupling 

The above splitting is further enhanced by the direct effect of the scalar isovector coupling, see Eg. I ll 11 1. Then the 
Schrodinger masses are 

ffi 

m* s (n,p) = (m* D ±f s p S3 ) — — — (18) 

m + (L s + L ) S y m =F (JppBl ~ fsPs3) 

The new term in the denominator will further contribute to the m* s (n) < m%(p) splitting since we must have 
fp > fs m order to get a correct symmetry parameter 04. The effect is larger at higher baryon densities because 

of the decrease of the scalar D53, due to the faster reduction of ps„- 

Actually in a non-relativistic limit we can approximate in Eq. (ll4t directly the energy e with the asymptotic kinetic 
energy leading to the much simpler relation: 

117 

-*- e- * m-Zo = >n* D - Ps+Zo)- (19) 



1 

In the case of asymmetric matter this leads to the more transparent relations 

• Only p meson coupling 

m* S ( n >P) = m *D - (?<s+'£<o)sym±fppB3 (20) 

• p + 8 coupling 

m* s (n,p) = m* Dsym - (L s + E ). (Vm ± (f p p B 3 - fsPs3) (21) 

In conclusion any RMF model will predict the definite isospin splitting of the nucleon effective masses ms(n)* < 
m s{p)* in me non-relativistic limit. We expect an increase of the difference of the neutron/proton mean field at high 
momenta, with important dynamical contributions that will enhance the transport effects of the symmetry energy. 



Beyond RMF: k-dependence of the Self-Energies 

Correlations beyond the Mean Field picture will lead to a further density and, most important, momentum depen- 
dence of the effective couplings and consequently of the Self-Energies. This can be clearly seen from microscopic 
Dirac-Brueckner (DBHF) calculations, see J ^flflfllll and refs. therein. Even the basic Fock correlations are implying 
a density dependence of the couplings, see 1131 . In fact from the very first relativistic transport calculations a momen- 
tum dependence of the vector fields was required in order to reduce the too large repulsion of the otpical potential at 
high momenta, see |28]. From the above discussion we can expect that, when the 5-meson is included 1 29], the sign of 
the D/rac-mass splitting will not be modified,i.e. m* Dn < m* Dp . In fact this is the result of the microscopic DBHF calcu- 
lations in asymmetric matter of refs.l9L lKXlllll . The problem is however that in correspondence the Schrodin ger mass 
splitting can have any sign. Indeed this is the case of the very recent DBHF analysis of the Tubingen group 1 1 ll 13(11 . 
For the 5-masses they get an opposite behavior m* Sn > m* s , with large fluctuations around the Fermi momenta, due to 
the opening of the phase space for intermediate inelastic channels. These results are actually dependent on the way the 
self-energies (mean fields) are derived from the full correlated theory and different predictions can be obtained, see 
the recent refs. 13 lLl32ll33ll . 

All that clearly shows that experimental observables are in any case strongly needed. 
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FIGURE 3. Energy dependence of the Dirac-Lane potential in the RMF picture (solid: NLp; dashed: NLpS) and in the 
phenomenologic Dirac Optical Model of Madland et al. jj^ljl/ , see text. 



The Dirac-Lane Potential 



In the non-relativistic limit of the Eq. (ll4> we can easily extract the neutron-proton mean optical potential using the 
isospin dependence of the self-energies 



^■O.q — ^0,symTfpPB3 , 
Z>s,q = ^s,sym ±/sPs3 , (22) 



With some algebra we get a compact form of the Dirac — Lane Potential 
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It is interesting to compare with the related discussion presented for the non-relativistic effective forces, mainly of 
Skyrme-like form. First of all we predict a definite positive E — slope given by the quantity f p which is actually not 
large within the simple RMF picture described here. This is a obvious consequence of the fact that the "relativistic" 
mass splitting is always in the direction ms(n)* < tns(p)*. The realistic magnitude of the effect could be different if we 
take into account that some explicit momentum dependence should be included in the scalar and vector self energies, 
1 28], as also discussed before. An interesting point in this direction comes from the phenomenological Dirac Optical 
Potential (Madland — potential) constructed in the refs. I37ll3^1 . fitting simultaneously proton and neutron (mostly 
total cross sections) data for collisions with a wide range of nuclei at energies up to 100 MeV. Recently this Dirac 
optical potential has been proven to reproduce very well the new neutron scattering data on 20& Pb at 96 MeV |39] 
measured at the Svendberg Laboratory in Uppsala. 

The phenomenological Madland — potential has different implicit momentum dependences (exp /log) in the self- 
energies, l37ll38ll . In Fig[3]we show the corresponding energy dependences for the Dirac — Lane potentials, compared 
to our NLp and NLp 8 estimations. When we add the 5-field we have a larger slope since the p — coupling should be 
increased. The slopes of the phenomenological potentials are systematically larger, interestingly similar to the ones of 
the Skyrme — Lyon forces. We note that in the Madland — potential the Coulomb interaction is included, i.e. an extra 
repulsive vector contribution for the protons. The first, not energy dependent, term of Ea. J23> is directly related to the 
symmetry energy at saturation, exactly like in the non-relativistic case, see Eq.Q. 

The conclusion is that a good systematic measurement of the Lane potential in a wide range of energies, and 
particularly around/above 100 MeV, would answer many fundamental questions in isospin physics. 



REACTION OBSERVABLES 



In order to directly test the influence of the S-effective mass splitting we will present results from reaction dynamics at 
intermediate energies analysed in a non-relativistic transport approach, of BNV type, see |5]. The I so —MD effective 
interaction is derived via an asymmetric extension of the GBD force, l34ll35ll . constructed in order to have a strict 
correspondence to the Skyrme forces, see also |36]. 
The energy density can be parametrized as follows: 

e = e + e kin + e(A',A") + e(B',B") + e(C\C") (24) 

where is the usual kinetic energy density and 

n 2 

e(A',A") = (A'+A"l5 2 )!— 
Po 

e(B'.B") = (B'+B"l3 2 ) ^j'p 
e(C,C") = C\J NN + J PP ) + C"J NP (25) 

Here ^ zx i are integrals of the form: 

J TT i= J dpdp'f T {r,p)f T ,{r,p')g(p,p') 

with g(p,p') = g(p — p') 2 This choice of the function g(p 7 p') corresponds to a Skyrme-like behaviour and it is 
suitable for BNV simulations. In this frame we can easily adjust the parameters in order to have the same density 
dependence of the symmetry energy but with two opposite n/p effective mass splittings. So we can separately study 
the correspondent dynamical effects, l36ll . In Fig. 0]we show the ^-dependence of the Lane potentials with the 
parametrizations corresponding to the two choices of the sign of the n/p mass splitting (shown in the insert for the 
/ = 0.2 asymmetry). The value of the splitting is exactly the same, only the sign is opposite. 
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FIGURE 4. Lane potential for the used parametrizations; small panel: the effective mass splitting, related to the slope of the 
Lane potential. 

The upper curve well reproduces the Skyrme-Lyon (in particular SLy4,Sly7) results. The lower (dashed) the 
SIH,SKM* ones. We clearly see the crossing at low energy, as expected from the previous discussion. 

From the figure we can immediately derive the expectation of very different symmetry effects for nucleons around 
lOOMeV kinetic energy, enhancement (with larger neutron repulsion) in the m* < m* p case vs. a disappearing (and even 
larger proton repulsion) in the m* > to* choice. So it seems naturai to look at observables where neutron/proton mean 
fields at high momentum are playing an important role. We will show results for fast nucleon emissions and collective 
flows at intermediate energies, in particular for high transverse momentum selections. 



Isotopic content of the fast nucleon emission 



We have performed realistic "ab initio" simulations of collisions at intermediate energies of n-rich systems, in 
particular 132 Sn+ 124 S« at 50 and 100, AMeV, central selection. Performing a local low density selection of the 
test particles (p < po /8) we can follow the time evolution of nucleon emissions (gas phase) and the corresponding 
asymmetry. In Fig. [5]we show the evolution of the gas asymmetry (7(f) = (N — Z) /A) for the 50AMeV reaction (the 
solid line gives the initial asymmetry). For both choices of the density stiffness of the symmetry term we clearly see 
the effects of the mass splitting, resulting in a reduced fast neutron emission when m* > m*. 
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FIGURE 5. Isospin gas content as a function of time in a central collision 132 Sn + 124 Sn at 50 AM eV for two opposite choices of 
mass splitting. Left: asy-soft EOS; right: asy-stiff EOS. 

In order to better disentangle the mass splitting effect and to select the corresponding observables, in Figs.|5J0we 
report the N/Z of the "gas" at two different times, t — 60/w/c, end of the pre-equilibrium emission, and t = IQOfm/c 
roughly freeze out time. We have followed the transverse momentum dependence, for a fixed central rapidity y^O) 
(normalized to projectile rapidity). In Fig.|6]we show the results without the mass splitting effect, i.e. only taking into 
account the different repulsion of the symmetry term (the intial average asymmetry is N/Z = 1 .56). 





2 


1 


.8 


N 




\1 


.6 


Z 






,4 




.2 



t=60fm/c2 



5 * 



Hi 



N 



6 

H T51.4 



200 400 
p, (MeV/c) 



• asy — soft 

• asy — siiff 

t= 100 fm/c 

.{*, 

: H 



200 400 
p, (MeV/c) 



FIGURE 6. Transverse momentum dependence of the neutron/proton ratio in the rapidity range \y^ \ < 0.3 for a central reaction 
I32 5n + 124 Sn at 50AMeV. Comparison between two choices of symmetry energy stiffness at different times. 

At early times (left panel of Fig. [6), when the emission from high density regions is dominant, we see a difference 
due to the larger neutron repulsion of the the asy-stiff choice. The effect is reduced at higher /?fS due to the overall 
repulsion of the isoscalar momentum dependence. Finally at freeze out (right panel) the difference is almost disap- 



pearing even for the more efficient Isospin Distillation of the asy-soft choice during the expansion phase, |5J] and refs. 
therein. 

When we introduce the mass splitting the difference in the isotopic content of the gas, for the larger transverse 
momenta, is very evident at all times, in particular at the freeze-out of experimental interest, see Fig. [7] As expected 
the m* > m* sharply reduces the neutron emission at high p t s. 
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FIGURE 7. Same as Fig.\6l for two opposite choices of mass splitting. 

We have repeated the analysis at higher energy, IQQAMeV, and the effect appears nicely enhanced, see Fig. [8] 
Similar results have been obtained in ref. 1 40] at 4-OOAMeV in a different I so — MD model, restricted to the 

m* > m* choice. In conclusion it appears that the isospin content of pre-equilibrium emitted nucleons at high transverse 

momentum can solve the mass-splitting puzzle. 



Differential collective flows and mass splittings 

Collective flows are expected to be very sensitive to the momentum dependence of the mean field, see 0| and 
refs. therein. We have then tested the isovector part of the momentum dependence just evaluating the Differential 
transverse and elliptic flows 

vl n - p \y,p t )^Vt 2 (y, Pt )-Vl 2 {y,p t ) 

at various rapidities and transverse momenta in semicentral (b/b max = 0.5) 197 Au + 197 Au collisons at250AMeV, where 
some proton data are existing from the FOPI collaboration at GSI l l4lll42il . 



Transverse flows 

For the differential transverse flows, see Fig. [5] the mass splitting effect is evident at all rapidities, and nicely 
increasing at larger rapidities and transverse momenta, with more neutron flow when m* < m* p 

Just to show that our simulations give realistic results we compare in lower right panel of Fig.|5]with the proton data 
of the FOPI collaboration for similar selections of impact parameters rapidities and transverse momenta. 
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FIGURE 8. Same as Fig. but for an incident energy of lOOAMeV . 



The agreement is quite satisfactory. We see a slightly reduced proton flow at high transverse momenta in the m* < m* 
choice, but the effect is too small to be seen fron the data. Our suggestion of the differential flows looks much more 
promising. 



Elliptic flows 

The same analysis has been performed for the differential elliptic flows, see Fig. [TO] Again the mass splitting effects 
are more evident for higher rapidity and tranverse momentum selections. In particular the differential elliptic flow 
becomes sistematically negative when m* < m*, revealing a faster neutron emission and so more neutron squeeze out 
(more spectator shadowing). 

In the lower right panel we also show a comparison with recent proton data from the FOPI collaboration. The 
agreement is still satisfactory. As expected the proton flow is more negative (more proton squeeeze out) when m* > m* . 
It is however difficult to draw definite conclusions only from proton data. 

Again the measurement of differential flows appears essential. This could be in fact an experimental problem due 
to the difficulties in measuring neutrons. Our suggestion is to measure the difference between light isobar flows, like 
triton vs. 3 He and so on. We expect to clearly see the effective mass splitting effects, maybe even enhanced due to 
larger overall flows, see o. 14311 . 
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